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Abstract 



We discuss the possible existence of new generic vacuum solutions of Robertson- 
Walker form in higher derivative gravity theories in four dimensions. These 
solutions illustrate how a dynamical coupling between very low and very high 
frequency modes can occur when the cosmological constant is small. 

The question of whether ultraviolet and infrared physics could in some way be 
coupled is intriguing. It is apparent that nonlinearities in some theory would have to 
play a role. But if we take the ultraviolet scale to be the Planck scale then Einstein's 
theory, although nonlinear, is inadequate to approach this question. We should at 
least consider a derivative expansion of some underlying theory of gravity, in which 
the Einstein action is simply the two derivative term. We advocate here that a 
study of generic solutions (no fine tuning of initial conditions) to generic actions (no 
fine tuning of coefficients in the derivative expansion) can shed some light on the 
possibilities, despite the purely classical nature of the approach. We find evidence 
for gravitational solutions with two widely disparate time scales. The period of a 
long time-scale oscillation is inversely related to the amplitude of a Planck time-scale 
oscillation, and both of these quantities are determined by the size of the cosmological 
constant. 

We consider a general gravitational action along with a free massless scalar field.[] 



We assume that all coefficients of the higher order terms are of order one in units of 
the Planck mass Mpi; e.g. a, b, c are of order Mp^. The only small dimensionless 
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^We shall see below that higher derivative terms involving the scalar are not important. 
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parameter is taken to be |A|/M|i. We consider the spatially flat metric, 

ds"^ = -dt^ + e^^^^ij dx'dx\ (2) 

along with (p = (j){t). The two gravitational field equations together imply the scalar 
field equation, + \B(f) = 0. 

An equation for B{t) alone can be obtained. Each term in this equation, besides 
the cosmological constant term, will have an even number of derivatives acting on 
some number of B{t)'s, where each B{t) has at least one derivative. The nonlinearity 
in the equation is essential for the solutions we study, but for \A\/Mpy <^ 1 it will 
turn out that all derivatives of B{t) are of order \f\A\ times the appropriate power 
of Mpi. Thus our solutions are hardly affected if we drop terms higher order than 
quadratic in B{t) in the equations. 

In addition we shall begin by dropping terms higher order than in (^. This 
truncation is not in any way justified, and in fact a problem will arise from this 
truncation. But the truncated theory serves to illustrate the basic dynamics, and 
we shall later see how the effects of higher order terms can eliminate the problem. 
The truncated equation can be put in the following dimensionless form involving 

m = m, 

A"" + A" + ciA'^ - c^A'A'" + c^A"^ = r] e\ (3) 

where ^ = t/ yjQa + 2h + 2c, = 4(3a + b + c)|A|, and r] = sign(A). In addition the 
gravity- scalar action fixes the nonlinear terms, 

ci = 3/2 , C2 = -9/2 , C3 = 3/2. (4) 

We have assumed that 3a + 6 + c > 0; this is the only combination of coefficients 
that survives for a conformally fiat metric, and we have absorbed it into the definition 
of our dimensionless expansion parameter e. The massless scalar contribution to the 
energy-momentum tensor corresponds to an equation of state p = wp with w = 1; for 
a more general contribution we find 

ci = 3(1 + w)/4: , C2 = -3 - 3w/2 , cg = 9/4 - 3w/4. (5) 

We shall examine solutions of (||) and treat (Q) and (||) as special cases. 

The de Sitter expansion A{^) oc ,^ is clearly a solution to @ when sign(ci) = 
sign(A). It can also be a generic solution in the sense that for a range of initial 
conditions the solution will tend towards the exponential expansion asymptotically. 
By a numerical analysis this is found to occur when 

- sign(2ci + C2) = sign(ci) = sign(A). (6) 
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The gravity- scalar action satisfies these conditions when A > 0, as does (|^) with 
w > —1. The initial conditions may be such as to produce fast oscillations initially, 
but these oscillations gradually die away as illustrated by the top curve in Fig. (1A).0 
Fig. (1) also illustrates a very different set of generic nonsingular solutions that 
exist when 

sign(2ci + C2) = sign(ci + C2 + C3) = sign(A). (7) 

By generic we again mean that if A'{0), A"{0), and A"'{0) are each in a respective 
range of values, all of order e, then one of these solutions results.0 We see that the 
gravity- scalar action satisfies these conditions as well, this time when A < 0. The 
same is true for a more general equation of state as long as w > 0. 

These new solutions are characterized by two scales of periodicity. One is an 
oscillation with a frequency of order unity in Planck units, as determined by the linear 
terms in (^), and an amplitude of order e. The other is a long time-scale oscillation 
with a frequency of order e and an amplitude of order unity. The amplitude of the 
fast oscillation is modulated with this same period. The complete solution is not 
strictly periodic since the phase of the fast oscillation in general changes after one 
period of the slow oscillation. 

The two scales of periodicity can be separated arbitrarily far since e can be arbi- 
trarily small. For example if A is a major component driving the present evolution 
of the universe, then the period of the slow oscillation can be of order the age of the 
universe. Given that the solution also requires a fast oscillation of order the Planck 
time, it is appropriate to refer to these solutions as UV-IR coupled. 

Two observations: 1) For a fine tuning of the initial conditions the amplitude of the 
slow oscillation can be brought to zero, leaving the fast oscillations only; these are the 
solutions studied in section I of 2) When e is not very small then the frequency of 
the two oscillations become similar, and some tuning of the initial conditions becomes 
necessary. 

Thus far we have considered an equation with up to 4-derivatives only. As we 
have alluded to above there is a problem with the UV-IR coupled solutions in this 
truncated theory. The other gravitational field equation determines an oscillating 
0^ that automatically satisfies the scalar field equation, but it also implies that 0^ is 
negative definite.^ This result is connected with having A < 0, since the contributions 

^The unphysical case in (j^) with A < would have exponential contraction as the generic solution. 

■^Derivatives of A(^) being of order e is related to the constants Ci being of order one. 

^In the de Sitter expansion example of Fig. (lA) the (p^ is also negative while asymptotically 
vanishing, but there is a different range of initial conditions that can produce positive (j)'^. A change 
in initial conditions cannot produce positive 0^ for the UV-IR coupled solutions to (||). 
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Figure 1: Solutions to equation (|^) with the coefficients in (^. In Figure A the top 
and bottom curves correspond to rje"^ = 10^'^/3 and —10^^/3 respectively. Figure 
B shows the continuation of the bottom curve in Figure A. Initial conditions have 
derivatives vanishing except for A"{0) = 0.035. Solutions with properties similar to 
these occur for a range of the initial conditions. 
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to the pressure p from A and the average 0^/4 tend to canceL Since 0^ contributes to 
p + p (unhke A) this problem with the truncated theory could be avoided if there was 
some other negative contribution to p+p. Higher order derivative terms involving the 
scalar field cannot help in this way since derivatives of 0^ are of order in Planck 
units (as can be seen from the scalar field equation), and 0^ itself is of order e"^. 

On the other hand all derivatives of A{C,) are of order e and thus contributions 
from higher powers of R can completely change the picture. They will add new linear 
and nonlinear terms to of the same order in e, respectively, as the ones in (^). We 
have studied the following equation with six derivative terms, where in general there 
is one additional linear term and three more quadratic terms. 

d^A""" + A"" + A" + -A'^ + -A' A'" + -A"^ 
i T 2 2 2 

+ d2A"'^ -d^A'A"" + d^A'A""' = r]e^ (8) 

For the two and four derivative terms we have inserted the values in (^. We find 
that for a range of coefficients of the new terms that UV-IR coupled solutions are 
once again generic, even when A > 0. 

The linear terms must support the fast oscillations; that is, the equation for the 
frequency obtained by inserting a sinusoid into the linear terms must have real roots 
only. Here it implies that < di < 1/4. Then to ensure that A and 0^ are positive, 
the quadratic coefficients must be such that ^2 + d^ is sufficiently large and positive 
and d4 is in some finite range that includes zero. It seems likely that these conditions 
can be satisfied for some range of the many coefficients that appear in the action 
with terms. But whatever the answer to that question, the question will have 
to be repeated for yet higher order terms. The answer can also be sensitive to the 
equation of state of other contributions to the energy momentum tensor. In the end 
a theory will either support the UV-IR coupled solutions or it will not; fine tuning is 
not necessary. 

Two examples of solutions to (^) with A > are displayed in Fig. (2), with 
the only difference being in the size of the initial conditions. For the larger initial 
derivatives. Fig. (2B) shows that the rapid oscillations can become irregular. This is 
presumably related to the existence of multiple rapidly oscillating modes associated 
with the linear terms, unlike the case in (^. 

We can at least address the meaning of the condition on the linear order terms at 
arbitrarily high order. The condition that the associated polynomial in the frequency 
has only real roots is just the requirement that there are no exponentially evolving 
modes among small fiuctuations that retain the form of the metric in (H). When 
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Figure 2: Solutions to equation (§) with coefficients di = .1, d2 = .2, ^3 = .2, 
= and rje'^ = 10~^/3. The initial conditions generating A and B have vanishing 
derivatives except for A"{0) = 0.005 and 0.028 respectively. Solutions exist for a 
range of the initial conditions, including the case where all derivatives in the initial 
conditions vanish. 
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extended to all orders in the derivative expansion of some underlying theory, this is 
the statement that the underlying theory is stable within this family of metrics. 

For a theory that supports the UV-IR coupled solutions with A > 0, there will 
also be the de Sitter solution A{^) oc ^. But in this case the de Sitter solution is 
unstable. Numerical analysis shows that the fast oscillations will tend to grow from 
any perturbation of the various derivatives of A{C,) to approach the order e amplitude, 
and eventually the large scale expansion will turn into a large scale oscillation. In other 
words, unlike the situation in (||), the de Sitter solution is not the generic solution. 
Of course in a different theory it may be that the de Sitter solution is generic, rather 
than the UV-IR coupled ones. For example in the de Sitter solution is generic for 

sufficiently large and positive. 

There is the related question of a possible decay || of the fast oscillations to 
energetic particles, occurring through the gravitational coupling. We see that this 
would involve a competition between the decay and the tendency for the oscillations 
to grow. And since the metric is conformally flat the decay would require conformal 
symmetry breaking, for example in the form of a particle mass, implying suppressed 
couplings of the oscillations to matter .0 

A decay to a de Sitter state would also assume that that is the preferred state. 
This is by no means obvious. Given the role of nonlinearities in the equations, it 
seems more likely that a more generic state is preferred where the various derivatives 
of A{C,) are all of the same order. Thus there may well be some kind of decay to a 
ground state, but that ground state could be a particular UV-IR coupled solution. 

All of this says nothing about the problem of why the cosmological constant 

is so small. This question becomes more interesting in some context in which the 

cosmological constant is a derived quantity. This may occur in higher dimensional 

theories 0, and in particular in a singularity-free warped Kaluza-Klein picture [Q 

where the effective 4D cosmological constant is a function of the amount of warping. 

In fact there are solutions [IH where time dependent oscillations coexist with the 

warping of a fifth dimension; and in this case the equations of the R + E? action 

yield 0^ > 0. But solutions having an exponential expansion as well were not found. 

This is understandable in light of the present results which indicate that the generic 

solutions are oscillatory on large time scales as well as short time scales. To actually 

demonstrate UV-IR coupled solutions within the warped Kaluza-Klein framework 

would require solving nontrivially coupled high-order differential equations. 

^ In higher dimensions the oscillating metric need not be conformally flat, resulting in larger 
couplings [||. 
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Nevertheless, UV-IR coupled solutions may well exist in such a framework. If 
such a solution is preferred for a given effective 4D cosmo logical constant, then the 
oscillation and its coupling to normal matter may provide a mechanism for the decay 
of the effective 4D cosmological constant. Note that the metric would have to adjust 
on both large and small time scales, since the relaxation would proceed through a 
series of UV-IR coupled solutions as e decreased. 
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